A BEALE-KATO-MAJDA BLOW-UP CRITERION FOR THE 3-D 
COMPRESSIBLE NAVIER-STOKES EQUATIONS 
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YONGZHONG SUN, CHAO WANG, AND ZHIFEI ZHANG 

Abstract. We prove a blow-up criterion in terms of the upper bound of 
the density for the strong solution to the 3-D compressible Navier-Stokes 
equations. The initial vacuum is allowed. The main ingredient of the 
proof is a priori estimate for an important quantity under the assumption 
that the density is upper bounded, whose divergence can be viewed as 
the effective viscous flux. 



1. Introduction 

In this paper, we consider the isentropic compressible Navier-Stokes sys- 
tem in three dimensional space. The system reads 

dtP + divipu) = 0, in (0, T) x Q, 

dtipu) + di\(pu ®u) - Lu + Vp = 0, in (0, T) x Q, 

together with the initial-boundary conditions 

(1.2) (p(t, x), u{t, x))\t=Q = (poix), uo(x)), in Q, 

(1.3) u(t,x) = 0, oniO,T)xdQ. 

Here Q. is either or a bounded domain in R^, p and u are the density and 
velocity of the fluid respectively, p = ap^ with y > 1 is the pressure. The 
Lame operator L is defined by 

Lu = //Am + (// + A)V div u, 

with constant viscosity coefficients p and A satisfying 

(1.4) n>0, 3i + 2yu>0. 

In the absence of vacuum for the initial density, the local existence of 
strong solution as well as the global existence of strong solution and weak 
solution with the initial data close to an equilibrium state were well devel- 
oped, see [|2Tll22ll24l[T4l[T0l l6l and references therein. The global existence 
of weak solution for large initial data was first solved by P. L. Lions in Il20l 
for y > |. E. Feireisl, A. Novotny and H. Petzeltova [13] extended Lions's 
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result to the case of 7 > |. S. Jiang and P. Zhang |[l8l[T9l| proved the global 
existence of weak solution for any y > 1 for the spherically symmetric or 
axisymmetric initial data. However, the regularity and uniqueness of weak 
solutions are completely open even in the case of two dimensional space. 
The only known result is the work of Kazhikhov and Vaigant fTl], where 
they proved the global existence of strong solution for the system (11.11 ) in 
Q = under the assumption that is a constant and X= pP with > 3. On 
the other hand, when the initial density is compactly supported, Z. Xin 1281 
proved that smooth solution will blow up in finite time in the whole space. 

To proceed we introduce some notations for the standard homogeneous 
and inhomogeneous Sobolev spaces. 

D' '(Q) =^ {u e LUa) : \\V'uhr,a^ < 00}, 

D^q(Q.) =^ {u 6 L\Q.) : WVuiy^ci) <oo and u = on dQ), 

When the initial vaccuum is allowed, the local well-posedness and blow- 
up criterion for strong solutions to the compressible Navier-Stokes equa- 
tions were established in a series of papers Cl[8l|9l| by Cho, Choe and Kim. 
Here we write down one of those results. 

Theorem 1.1. Let f2 be a bounded smooth domain or K.^ and q 6 (3,6]. 
Suppose that po > and belongs to W^''^{Q) n H\Q) n L^iQ), Uq e D^^iQ) n 
D^(Q.) with the following compatibility condition satisfied, 

(1.5) - i^Auo -(ju + A)V div uq + Vpipo) = yfpog, 

for some vector field g 6 L^(Q.). Then there exist a time T e (0, 00] and a 
unique strong solution (p, u) to (|l.ll) such that 

p 6 C([0, T), n W^'\Q.)), u 6 C([0, T), 0^(0)) n L^iO, T; D^'\Q.)). 

Moreover, let T* be a maximal existence time of the solution. If T* < 00, 
then there holds 

(1.6) limsup(||p(0llwi.7(n) + l|M(Ollz)'(n)) = ^■ 

Since the initial vacuum is allowed, it is then important to investigate the 
possible blow-up mechanism of the solution. In their recent works [[T5l[T6l . 
X. Huang and Z. Xin established a Beale-Kato-Majda blow up criterion for 
the above strong solution. More precisely. 
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Theorem 1.2. Assume that the coefficients of the operator L satisfies (11.41) 
and moreover, 

(1.7) A<ltx. 

Let (p, u) be the strong solution constructed in Theorem 11.11 and T* be a 
maximal existence time. If T* < oo, then 

(1.8) lim ||VM||z,i(o,r,L"(n)) = ^o- 

Recently, J. Fan, S. Jiang and Y. Ou fV2\ also obtained a similar result for 
the compressible heat-conductive flows. On the other hand, for the 2D com- 
pressible Navier-Stokes equations in T^, B. Desjardins IfTTI proved more 
regularity of weak solution under the assumption that the density is upper 
bounded; Very recently, L. Jiang and Y. Wang [17], Y. Sun and Z. Zhang 
||26ll obtained a blow-up criterion in terms of the upper bound of the den- 
sity for the strong solution. In [|26ll . the initial vacuum is allowed and the 
domain includes the bounded domain. Note that the L'(0, T, L°°(Q)) bound 
for Vm immediately implies the upper bound for the density p. 

The purpose of this paper is to obtain a Beale-Kato-Majda blow-up cri- 
terion in terms of the upper bound of the density for the 3-D compressible 
Navier-Stokes equations. Our main result is stated as follows. 

Theorem 1.3. Assume that (p, u) is the strong solution constructed in The- 
orem ll.ll Let yu, A be as in Theorem I 1.21 and T* be a maximal existence time 
of the solution. If T* < oo, then we have 

(1 .9) lim sup ||p(OllL~(o,r;t"(n)) = ^o- 

Remark 1.4. This result seems surprise, if we compare with the incom- 
pressible Navier-Stokes equations where the density is a constant. It is 
well-known that if we have some kind of control for the pressure, the Leray 
weak solution is in fact smooth for the incompressible Navier-Stokes equa- 
tions, see [T]. For the compressible Navier-Stokes equations, the pressure 
is determined by the density, the bound of the density thus implies a bound 
for the pressure. From this viewpoint, our result seems natural. 

Remark 1.5. In a forthcoming paper, we will extend similar result to the 
compressible heat-conductive flows. 

Let us conclude this section by introducing the main idea of our proof. 
First of all, if the density is upper bounded, we can obtain a high integra- 
bility of the velocity, see Lemma 13.11 This bound can be used to control 
the nonlinear term. The trouble is to control the density, which satisfies a 
transport equation. In order to propagate the regularity of the density, it is 
necessary to require that the velocity is bounded in L}{Q, T; W^'^'iQ.)). On 
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the Other hand, we have to obtain a priori bound of Vp in order to prove 
u 6 L'(0, T; W^'"'{Q.)). To overcome this difficulty, we introduce an im- 
portant quantity w defined by w = w - v, where v is the solution of Lame 
system 

//Av + (A+ iu)V div V = Vp(p) in Q, 
v(x) = on dCl. 

def 

In the case of Q = R^, (i + 2yu)divw = (A + 2yu)divM - p = G. It is 
well known that G is called the efl'ective viscous flux, which plays an im- 
portant role in the existence theory of weak solution. A key point is that 
we can obtain the better regularity of w than u under the only assump- 
tion that the density is upper bounded. More precisely, we proved that 
V^w e L^(0, T\ L^{Q.)), which combined with the bound of the density im- 
pUes that Vu e L^(0, T;L"'(Q) + L'^(0, T; BMO(Q)). This bound stiU does 
not imply that Vu is bounded in L^(0, T ; L'^{Q.)). We need to introduce the 
second key ingredient: a logarithmic estimate for Lame system. Then the 
result can be deduced by combining the above two estimates into the energy 
estimates for the density. 



2. Preliminaries 

Consider the following boundary value problem for the Lame operator L 
fiAU + (m + A)V div U = F, in a, 
U{x) = 0, on dn. 

Here U = {U\, U2, C/3), F = (Fi, F2, F3). It is well known that under the 



(2.1) 



assumption (II. 4|) . (12.11) is a strongly elliptic system. If F 6 W~^'^{Q.), then 
there exists an unique weak solution U 6 Z)q(Q). We begin with recalling 
various estimates for this system in spaces. 

Proposition 2.1. Let q e (1, 00) and t/ be a solution of (12.11) . There exists a 
constant C depending only on A, fx, q and Q. such that the following estimates 
hold. 

F e Uia), then 



(2.2) 



iy2.9(n) ^ if^^is abounded domain. 

(2)If F 6 F = div / with / = (/,)3x3,i;7 6 L^Q)), then 



(2 3) I - 

\ ||i7|liyi.9(n) ^ C'||/||L9(n); if is a bounded domain. 
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(3) If F = div / with = dkh]. and h]. e w!^'\a) for i, j,k= 1,2, 3, then 
(2.4) WUWma) < C||/i||L,(n). 

Proof. In the case when 11 is a bounded domain, the estimates (|2.2I) and 
(12.31) are classical for strongly elliptic systems, see for example [3J. The 
estimate (12.41) can be proved by a duality argument with the help of (|2.2I) . 
In the case of Q = R^, one can give an explicit representation formula for 
the solution as follows. Taking divergence on both sides of (12.11) . one finds 

diwU = — i — A'MivF. 
A + 2fi 

Substituting this into (12.11) gives us 

AU = -F VA"^ divF. 

yU + 2/i) 

Denote the Riesz transform = (i?i,i?2,^3) = VA"'''^. Then 

1 A + Li 

AU = -F —R{R ■ F). 

fi + 2//) 

Hence for /, j,k = 1,2,3, 

dijUk = -RiRiFk RiRiRkiR ■ F). 

The classical L^(R^)-boundedness for Riesz transform gives 
\\D'U\\u,^w^ < C(^)4T^lli^llL.(R3). 

IJ.{A + 2yU) 

Similar argument gives the estimates (12.31) and (12.41) . □ 
We need an endpoint estimate for L in the case q = oo. Let BMO(Q.) 
stand for the John-Nirenberg's space of bounded mean oscillation whose 
norm is defined by 



with 



def 

WfllBMO^a) = UIIlHO.) + [f]BMO, 



[fhMO(n)^= sup f \f(y) - fii,(x)\dy, 

xeSl,re{0,d) J 

fa,x, = j f(y)dy = ^ J f(y)dy. 



n,(.r) n,.{x) 

Here Q.r(x) = Br(x) n Q, Br{x) is the ball with center x and radius r and d 
is the diameter of Q.. \Q.r(x)\ denotes the Lebesque measure of Q.r{x). Note 
that 

[f]BMO(a) ^ 2||/||ioo(j^). 
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Proposition 2.2. If F = div/ with / = (fijhxsjij 6 L°°(Q) n L^{Q), then 
VU e BMO(Q.) and there exists a constant C depending only on A,^ and 
such that 

(2.5) l|Vt/||fiMO(n) < C (||/||L»(n) + UhHii)) • 

Proof. When Q is a bounded domain,the estimate (12.51 ) can be found in 
|[T1 for a more general setting. Now if ^1 = R.^ we use the representation 
formula for Vi7. Since 

At/ = - div / — ^ VA-i div div / = - div / — ^ VG, 

/i fi{A + 2yU) yU yU(/l + 2/J.) 

with G = Z-;=i ^/^y/i- For = 1, 2, 3, 

IX iu(A + 2/j.) 

By the Fefferman-Stein's classical result on fiMO-boundedness of singular 
integral operators Il251 . there exists an absolute constant C > such that 

2/1 ~\~ ^tju 

[Vt/]BMO(R3) ^ C— — T— -||/||L»(n)- 
IX(A + 2yU) 

This inequality combined with (|2.3I) with ^ = 2 yields (|2.5I) . □ 
In the next lemma, we will give a variant of the Brezis-Waigner's inequal- 
ity [l5j|. To our knowledge, such a kind of inequality was first established 
in [|23]| in the case of f2 = R^. For the reader's convenience, we will give a 
proof in the case when Q is a bounded Lipschitz domain, see also [i26l . 



Lemma 2.3. Let Q = 1.^ or be a bounded Lipschitz domain and / e 
W^'''{Q.) with q e (3,oo). There exists a constant C depending on q and 
the Lipshitz property of Q. such that 

(2.6) Wfh-in) < C(l + ||/||BMO(n)ln(e + WVfWma))) ■ 

Proof. First note that for a Lipschitz domain, the following so-called A- 
property holds: 

There exist two constants A > 1 and ro e (0, d) such that for any r 6 
(0, ro) and x e Q., 

Without loss of generality we assume rg < 1. 

First we give an estimate for |/n,(,:)| with < r < ro and x 6 Q. If r > ^ro, 
then 
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If r < jro, then there exists some integer k> \ such that 
^<r<'f^,k<C(l+\lnr\). 

Denoting Qj = Q.2jr(x) for j = 0, 1, • • ■ , fc, we have 

k 

|/n,(A-)| < ^ |/n,_, - /n, | + |/nj 

< yf |/(y)-/n,ky + C||/|^„) 



< 2^Ay f |/0;)-/^Jj3; + C 



LHa) 



< Ck[fhMOin) + C\\f ly^n) < C(l + |ln r\)\\f\\sMO(,i). 
We conclude that there exists a constant C = C(A, ro, AO such that 

|/aw|<C(l+|lnr|)||/||BMO(n), 

which together with Sobolev embedding theorem in a Lipschitz domain [|2l 
ensures that for any fixed x e CI and small enough e > we have 

\f{x)\ < \f{x) - + \fnAx)\ < C(e'-f 11/11 wi..(n) + (1 + I lne|)||/|| 

A suitable choice of s yields the inequality (12.61) . □ 
In the subsequent context we will use L~'F to denote the unique solution 



U of the Lame system (12.11) . 

3. A PRIORI ESTIMATES FOR THE EFFECTIVE VISCOUS FLUX 

In what follows, we assume that (p, u) is a strong solution of (11.11) in 
[0, T) with the regularity stated in Theorem 1 1.1[ 

Standard energy estimates yields that for any t e [0, T), 



\l\q.) ^ IIPollLHn)^ 
llp(OIII.(n) + llpl"l'(OllL'(n) + IIVw||^2((o,,)xn) 
< Cdlpoll Lr(n) + IIPolwoPllLHn))- 
Note that by the assumption on p^, uq, 

llPoillr(n) < IIPollI^Jn)lij^ollLi(n)' IIPol"oPllL'(n) < WpoWtyHiDlMllef^^y 
We thus have the following bounds 

(3.1) llpllL"(0,7';Li(n)), II \^M|lL~(0,r;L2(n)), l|VM||L2(0,r;L2(n)) < C. 

Here C depends only on/u,A,y,a and po, mo- 
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In what follows the dependence of the constant C on jx, X, y, a and Q. will 
not be mentioned. 

The following lemma is the first key step, whose argument comes from 
m and [fT6ll. 

Lemma 3.1. Assume that n <1A and the density p satisfies 

(3.2) llpllL~(o,r;L»(n)) < M. 

There exists r 6 (3, 6) such that p\uY 6 L~(0, T;L\a)) with 

Here C depends on T, ||pollL~(n), IIVi/ollL^cn), 

Proof. Multiplying the second equation of (|l.ll) by r\uY''^u, and integrating 
the resulting equation on Q to obtain 

a. £1 
+r(r - 2) [/x\ur^\V\uf + {A + fi){div uM'-^ ■ V\u\) dx 



(3.3) = ^ rp(p)di\(\u\'' ^u)dx 

By using the fact |Vm| > IV|mII, the term in the second integrand can be 
estimated from below by 

r\ur^[lx\Vu\^ + {A+ iS){d\\uf + (r - 2)n\V\u\\^ 

-(i+yU)(r-2)|VN||diVM|] 

> r\ur^\pi\VuY + {A + div u - '-^\V\u\\f -(A + ^)^J-:^\V\uf 

+(r-2)//|V|M|p] 

> r\ur^[ju\Vu\^ + (r - 2){m - (A + fi/-^)\V\u\Y] 

Recalling that A < Iji, there exists r e (3, 6) such that the last term is greater 
than 

c\u\'-^\Vu\^. 

On the other hand, because of ||pIIl~ < M, we find that the right-hand side 
of (|3.3I) is controlled by 

C ^ p'^\ur^\Vu\dx < e J \ur^\Vu\^dx + ^ J pHdx 



n. n 
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Taking e = | to yield that 



'■-2 



^ J" p\u\''dx < C J p\u\^dx 



which together with the following bound 

IIPolMoniLi(n) < IIPoll^5i,(^jllMolli6(n) ^ C'||poI|^^^^^||VmoII22(j^), 

implies the desired estimate. □ 

def 

Now for each t 6 [0, T), we denote v{t,x) = L'^Vpip). That is, v{t) is 
the solution of 

'//Av + (/I + yu)Vdivv = Vp(p) in Q., 
v(t, x) = on da. 



(3.4) 



Thanks to Proposition 12. 1[ for any q e (l,oo), there exists a constant C 
independent of t such that 



(3-5) „^2 



l|Vv(OllL.(n) < C\\p(pm\Lmh 



Now let us introduce an important quantity 

W = U -V, 

whose divergence can be viewed as the effective viscous flux. 

An important observation is that this quantity possesses more regular- 
ity information than u does under the assumption that the density is upper 
bounded. More precisely, 

Proposition 3.2. Under the assumption (13.21) . we have 

(3.6) ||Vw||i~(o,r;L2(£i)), llp^5(W||i2((o,T)xn), IIV^w||i,2((o,7-)xn) < C. 

Here the constant C depends on ||pollL~(f2). IIVwollL^cn), M, T. 

Proof. By using the continuity equation, we find that w satisfies 

• pdfW - luAw - (i + yu) V div w = pF, in (0, T) x a, 
w(t, x) = Oon [0, r) X da, w(0, x) = Wq(x), in Q, 

with Wq(x) = Uo(x) + Vo{x) and 



(3.7) 



F = -u-Vu-L-'Vid,p(p)) 

= -u ■ Vu + L'^V div[p(p)u] + L~^V[(pp'(p) - pip)) div u]. 
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Multiplying the first equation of (|3.7I) by dtW and integrating the resulting 
equation over Q to obtain , 

//|Vvvp + (/l + yu)|divwpJjc + J p\dtw\^dx = J pF ■ dtwdx, 
n n n 

which together with Holder inequality and Young's inequality gives 

n n 

(3.8) ^^ll^^^ll^.(n)- 

Now let us estimate || y/pFW^^if^^y We get by Lemma [XT] and (13.51) that 

< c||p^.|L,.)(||Vw||^ . + 

< C,||Vw||i2(n) + e||V2>v|L2(n) + C. 
Here we use the interpolation inequality 

ll/llL.(n) < CMhHn) + e||V/|L2(n), 2 < q < 6. 
We infer from Proposition 12.11 that 

II ^L-'Vdiv[p(p)M]||i2(n) < CIip(p)Mili2(n) < C|| yfpuiy^a) < C, 
II y/pL~^V(pp' - p)divM||i,2(n) 

< li^^lb(a)||L-iV(pp'-p)div 

< C\\VL-'V(pp' -p)diwuh2^a) < Cimy^ay 
Consequently, for 6 > to be determined later, 

(3.9) II ^^F||^,^^^ < fllV^wll^,^^^ + C,(l + ||Vw||i,(^^ + WVuWl,^^;). 
Noting that Lw = pdtW - pF, we get by using Proposition 12.11 again that 
W'^'^Hm ^ Ci\\pdM?^2,^, + WpFWl.,^,) < C(|| y/pdMll^,^, + II ^^pF\\l2,^,l 
which implies by taking e = in (13.91) that 

II ^pwIh^, ^ ^11 yfpdM\l2,a, + c(i + livwii^,^^^ + llVi^ii^,^^^). 

Substituting this estimate into (IX8l) and noting that l|VM(0ll^2(n) ^ 
the estimate (|3.6I) follows from Gronwall's inequality. □ 
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Corollary 3.3. Under the assumption (13.21) . we have 

\\^n\\L°'{OJ-LHn)), l|w|lL"(0,7-;L6(n)), \\^H\\LH0,T-Li(n)) ^ C, 

for any q e [2, 6]. 

Proof. This can be deduced from Proposition 13. 2[ (13. 5|) and Sobolev em- 



4. High order a priori estimates for the effective viscous flux 

In this section, we will give high order regularity estimates for w. This is 
possible if the initial data (po, "o) satisfies the compatibility condition (|1.5I) . 
We still assume that (p, u) is a strong solution of (11.11) in [0, T) and satisfies 
(13. 21) . The energy estimates in this section are motivated by the calculations 
of D. Hoff^lUl. 

We begin by introducing some notations. For a function or vector field(or 
even a 3 x 3 matrix) x), the material derivative / is defined by 



and div(/ ® u) = Yj]=\ djifuj). For two matrices A = (a,y)3x3 and B = 
(l>ij)sx3, we use the notation A : B = 2^y=i ^O'^/y ^i^d AB is as usual the 
multiplication of matrix. 

We rewrite the second equation of (11.11) as 



By taking the material derivative to the above equation and using the fact 
f = ft + div{fu) - f div u, we obtain 

piif + pu-Vu-\- Vp, + div(Vjc» (8) u) 
(4.1) = //[Awf + div(AM (8» u)] +(A + ju)[VdivM, + div((VdivM) ® u)]. 

Multiplying (14.11) by u and integrating on Q to obtain 



bedding theorem. 



□ 



pii + Vp(p) - Lu = 0. 




a 



n 



(4.2) 




-{A + fj.) I M • ((V div Ut) + div((V div u) (8) u)))dx 



a 
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The //-term can be calculated as follows. 

- J" M • (Am, + div(AM (8» u))dx = [Vu : Vut + m (S> Am : Vm] dx 
a n 

= J [|VmP - V(m ■ Vm) : Vm + m (8) Am : VmJJx 

n 

= J [|Vm|^ - {(VuVu) + (m • V)Vm) : Vm - V(m • Vm) : Vujdx 

n 

= J [|VmP - (VuVu) : Vm - div(VM ®u):Vu 
a 

-(VuVii) : Vm - ((m • V)Vm) : Vu\dx 

= J [|VmP - (VuVu) : Vm + ((m • V)Vm) : Vm 
a 

-(VuVu) : Vm - ((m • V)Vm) : VMjjx 

/r 



> J -IVmP-CIVm^ 

To estimate the (A + yu)-term of (14.21) . note that 

div((V div u)®u) = V(m • V div u) - div(div mV O m) + V(div uf, 
div M = div Mf + div(M • Vm) = div m, + m ■ V div m + Vm : (Vm)'. 

Here A' means the transpose of matrix A. We have 

- J" M ■ [v div M, + div((V div u) ® u)^dx 
n 

= ^ ^ div M div u, + div m(m • V div m) 

n 

- div uiyu)' : Vm + div M(div M)^j 

= J" [l div Mp - div mVm : (Vm)' - div m(Vm)' : Vm + div M(div M)^j^^x 
n 

> J [^IdivMp- ^|VmP-C|Vm|^]Jx 
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We continue to estimate the pressure term. 
Pt div ii + (u • Vm) ■ Vpdx 



I' 



n 

= p'ip)Pt div ii + (u ■ Vw) • Vpdx 

Si 

= J" -pp'ip) div M div M - (m • Vpip)) div u + (u ■ Vm) • Vpdx 
a 

= J" -pp'ip) div udivu + div((div u)u) - div{{u ■ Vu))^dx 
a 

= -pp'ip) div u div ii + div u div u - (Vm)' : V^jj^ 

n 

where we used the assumption (13.21) and Corollary 13.31 in the last two in- 
equalities. 

Substituting those estimates into (|4.2I) yields 

J p\u\^dx + p J \Vu\^dx + iA + p) J\di\u\^dx 
n n. n 

(4.3) <C J Wufdx + CWVuiy^ny 



a 



To conclude the estimate by Gronwall's inequality, we will use the term 
II y/pii\\L^(D.) to control ||VM||i4(f|). Thanks to the definition of w, we know that 
w satisfies 

(4.4) pAw + iA + p)V div w = pu in Q., 

with the zero boundary condition. We get by Proposition 2. 1 that 

l|V^w||L2(n) < C||pM||L2(n) < C|| ^M||z.2(n), 

which together with the interpolation inequality, Corollary 13. 3[ and Propo- 
sition 2.1 leads to 

||VM||^4(ji) ^ IIVM|lL2(n)IIVM||j6(jj) < C||VM||i6(n)||VM||2,^^^ 

< C||VM||J,(^/||Vw|L<,(n) + \\Vv\y(ci)) 

< C||Vm||^,(^^(i + II ^M||L2(n)). 
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Substituting this estimate into (l431)and noting that l|Va(0ll^6(j^) e ^HO, T) 
by Corollary 3.3, we get by Gronwall's inequality that 

(4.5) ^ p\u\^dx + J \Vu\^dxdt < C, 



n n 



with C depending only on T, M and po, Uo,g. Here we used the compatibil- 
ity condition (11.51) . 

With the help of Sobolev embedding theorem and using the equation (14.41) 
again, we deduce from (14.51) that 

Proposition 4.1. Under the assumption (3.2), we have for all 2 < q < 6, 
(4.6) l|V>v||i2(o,7-.i"(n)), W^w\\LHo,T-Lnm < C, 

with the constant C depending on q, M, T and po, Uo, g- 

5. Proof of Theorem 1 1.3 1 

Now we are in position to prove Theorem 11.31 We will prove it by the 
contradiction argument. Assume that T* < oo and 

sup ||p(5)||z.»(n) < oo. 

.ve[0,7-) 

By Theorem ll.il it suffices to show that 

(5.1) sup ||Vp(5)||L,(n) < oo. 

se[0,T') 

Taking the derivative with respect to x for the first equation of (11.11 ) to 
obtain 

(5.2) dtVp + (u ■ V)Vp + VwVp + div uVp + pV div u = 0. 
In the following estimates we will use 

(5.3) WVMmm ^ QWphnn), 

l|Vv|L~(f^) < c(l + \\Vv\\BMO(iiMe + \\V\\\mii))) 

< C(l + Wph-nLHci) He + l|Vp||L,(n))) 

(5.4) < c(l+ln(e + ||Vp|b(n))) 

with the second estimate followed from Proposition 12. 1[ 12.21 and Lemma 
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Multiplying (|5.2I) by <?|Vp|^ ^Vp and integrating the resulting equation on 
Q., we obtain 



j^J \Vp\"dx ^ cj IVuWVpl'idx + q J p\Vdwu\\Vpr^dx 
£1 ii a 



< C(||Vw|L»(n) + ||Vv|L»(n))||Vp|| 

+C{\\VM\m^) + W'v\\mn))\mV 



1 

Li(0.) 



from which and (I5.3I) - (I5.4I) . we infer that 



j^^WpVdx < C(l + ||Vv||™ + ||Vw|L»(n))||Vp|| 

72,, ,11 IIY7^ll?-l 



1 

Li{0.) 



+ciivviL,(n)iivpii^;;^) 

< C(l + ||Vw||w.,,(n) + ln(e + ||Vp|L,(n)))l|Vp||^,(^3 

+l|V'>v|L,(n)l|Vp||^4. 

Note that \\Vw\\^^U(n) 6 L^(0,T*) by Proposition 1411 Then by Gronwall's 
inequality, we conclude the proof of (15.11 ) and hence Theorem 1 1.31 □ 

Acknowledgement 

The authors thank Professors Song Jiang, Changxing Miao, Zhouping 
Xin and Ping Zhang for their profitable discussion and suggestions. This 
work was done while Yongzhong Sun and Zhifei Zhang were visiting the 
Momingside Center of Mathematics in CAS. We would like to thank the 
hospitality and support of the Center. Yongzhong Sun is supported by NSF 
of China under Grant 10771097. Zhifei Zhang is supported by NSF of 
China under Grant 10990013. 



References 

[1] P. Acquistapace, On BMO regularity for linear elliptic systems, Ann. Mat. Pura Appl., 
161 (1992), 231-269. 

[2] R. A. Adams, Sobolev spaces. Pure and Applied Mathematics, Vol. 65, Academic 

Press, New York-London, 1975. 
[3] S. Agmon, A. Douglis and L. Nirenberg, Estimates near the boundary for solutions 

of elliptic partial differential equations satisfying general boundary conditions. I. 

Comm. Pure Appl. Math. 12(1959), 623-727. 
[4] L. C. Berselli and G. P. Galdi, Regularity criterion involving the pressure for the 

weak solutions to the Navier-Stokes equations, Proc. Amer Math. Soc, 130(2002), 

3585-3595. 

[5] H. Brezis and S. Wainger, A note on limiting cases of Sobolev embeddings and con- 
volution inequalities. Comm. R D. E., 5(1980), 773-789. 



1 6 YONGZHONG SUN, CHAO WANG, AND ZHIFEI ZHANG 

[6] Q. Chen, C. Miao and Z. Zhang, Global well-posedness for the compressible Navier- 
Stokes equations with the highly oscillating initial velocity, to appear in Comm. Pure. 
Appl. Math. 

[7] Y. Cho, H. J. Choe and H. Kim, Unique solvability of the initial boundary value 

problems for compressible viscous fluids, J. Math. Pures Appl., 83(2004), 243-275. 
[8] Y. Cho and H. Kim, On classical solutions of the compressible Navier-Stokes equa- 
tions with nonnegative initial densities, Manus. Math. 120(2006), 91-129. 
[9] J. H. Choe and H. Kim, Strong solutions of the Navier-Stokes equations for isentropic 
compressible fluids, J. Diff. Equa., 190(2003), 504-523. 

[10] R. Danchin, Global existence in critical spaces for compressible Navier-Stokes equa- 
tions. Invent. Math., 141(2000), 579-614. 

[11] B. Desjardins, Regularity of weak solutions of the compressible isentropic Navier- 
Stokes equations. Comm. P D. E., 22(1997), 977-1008. 

[12] J. Fan, S. Jiang and Y. Ou, A blow-up criterion for the compressible viscous heat- 
conductive flows, Annales de I'lnstitut Henri Poncare-Analyse non lineaire, online. 

[13] E. Feireisl, A. Novotny and H. Petzeltova, On the global existence of globally defined 
weak solutions to the Navier-Stokes equations of isentropic compressible fluids, J. 
Math, fluid Mech., 3(2001), 358-392. 

[14] D. Hoff, Global solutions of the Navier-Stokes equations for multidimensional com- 
pressible flow with discontinuous initial data, J. Diflr. Equa., 120(1995), 215-254. 

[15] X. Huang and Z. Xin, A blow-up criterion for the compressible Navier-Stokes equa- 
tions, arXiv:0902.2606 

[16] X. Huang and Z. Xin, A blow-up criterion for classical solutions to the compressible 
Navier-Stokes equations, arXiv: 0903.3090. 

[17] L. Jiang and Y. Wang, On the blow up criterion to the 2-D compressible Navier-Stokes 
equations. Preprint, 2009. 

[18] S. Jiang and P. Zhang, Global spherically symmetric solutions of the compressible 
isentropic Navier-Stokes equations. Comm. Math. Phys., 215(2001), 559-581. 

[19] S. Jiang and P. Zhang, Axisymmetric solutions of the 3-D Navier-Stokes equations 
for compressible isentropic flows, J. Math. Pure Appl., 82(2003), 949-973. 

[20] P. L. Lions, Mathematical Topics in Fluid Mechanics, V.2, Compressible Models, 
Clarendon Press, Oxford, 1998. 

[21] A. Matsumura and T. Nishida, The initial boundary value problems for the equations 
of motion of compressible and heat-conductive fluids. Comm. Math. Phys., 89(1983) 
445-464. 

[22] J. Nash, Le probleme de Cauchy pour les equations differentielles d'un fluide general. 
Bun. Soc. Math. France, 90(1962), 487-497. 

[23] H. Kozono and Y. Taniuchi, Limiting case of the Sobolev inequality in BMO, with 
application to the Euler equations. Comm. Math. Phys., 214(2000), 191-200. 

[24] V. A. Solonnikov, The solvability of the initial-boundary value problem for the equa- 
tions of motion of a viscous compressible fluid, J. Sov. Math., 14(1980), 1 120-1 133. 

[25] E. M. Stein, Harmonic analysis: Real-variable methods, Orthogobality, and Oscilla- 
tory integrals, Princeton Univ. Press, Princeton, 1993. 

[26] Y. Sun and Z. Zhang, A blow-up criterion of strong solution for the 2-D compressible 
Navier-Stokes equations. Preprint, 2009. 

[27] V. A. Vaigant and A. V. Kazhikhov, On the existence of global solutions of two- 
dimensional Navier-Stokes equations of a compressible viscous fluid, Sibirsk. Mat. 
Zh., 36(1995), 1283-1316; translation in Siberian Math. J., 36(1995), 1108-1141. 



COMPRESSIBLE NAVIER-STOKES EQUATIONS 



17 



[28] Z. Xin, Blowup of smooth solutions to the compressible Navier-Stokes equation with 
compact density, Conmi. Pure Appl. Math., 51(1998), 229-240. 

Department of Mathematics, Nanjing University, 210093, P. R. China 
E-mail address: simyongzhong@163 . com 

Academy of Mathemahcs & Systems Science, CAS, Beijing 100190, P. R. China 
E-mail address: wangchao@amss .ac.cn 

School of Mathematical Sciences, Peking University, 100871, P. R. China 
E-mail address: zf zhangOmath . pku . edu . cn 



